An autoantibody signature predictive for multiple sclerosis.

Although B cells are implicated in multiple sclerosis (MS) pathophysiology, a predictive or diagnostic autoantibody remains elusive. In this study, the Department of Defense Serum Repository (DoDSR), a cohort of over 10 million individuals, was used to generate whole-proteome autoantibody profiles of hundreds of patients with MS (PwMS) years before and subsequently after MS onset. This analysis defines a unique cluster in approximately 10% of PwMS who share an autoantibody signature against a common motif that has similarity with many human pathogens. These patients exhibit antibody reactivity years before developing MS symptoms and have higher levels of serum neurofilament light (sNfL) compared to other PwMS. Furthermore, this profile is preserved over time, providing molecular evidence for an immunologically active preclinical period years before clinical onset. This autoantibody reactivity was validated in samples from a separate incident MS cohort in both cerebrospinal fluid and serum, where it is highly specific for patients eventually diagnosed with MS. This signature is a starting point for further immunological characterization of this MS patient subset and may be clinically useful as an antigen-specific biomarker for high-risk patients with clinically or radiologically isolated neuroinflammatory syndromes.

High-throughput complement component 4 genomic sequence analysis with C4Investigator.

The complement component 4 gene loci, composed of the C4A and C4B genes and located on chromosome 6, encodes for complement component 4 (C4) proteins, a key intermediate in the classical and lectin pathways of the complement system. The complement system is an important modulator of immune system activity and is also involved in the clearance of immune complexes and cellular debris. C4A and C4B gene loci exhibit copy number variation, with each composite gene varying between 0 and 5 copies per haplotype. C4A and C4B genes also vary in size depending on the presence of the human endogenous retrovirus (HERV) in intron 9, denoted by C4(L) for long-form and C4(S) for short-form, which affects expression and is found in both C4A and C4B. Additionally, human blood group antigens Rodgers and Chido are located on the C4 protein, with the Rodger epitope generally found on C4A protein, and the Chido epitope generally found on C4B protein. C4A and C4B copy number variation has been implicated in numerous autoimmune and pathogenic diseases. Despite the central role of C4 in immune function and regulation, high-throughput genomic sequence analysis of C4A and C4B variants has been impeded by the high degree of sequence similarity and complex genetic variation exhibited by these genes. To investigate C4 variation using genomic sequencing data, we have developed a novel bioinformatic pipeline for comprehensive, high-throughput characterization of human C4A and C4B sequences from short-read sequencing data, named C4Investigator. Using paired-end targeted or whole genome sequence data as input, C4Investigator determines the overall gene copy numbers, as well as C4A, C4B, C4(Rodger), C4(Ch), C4(L), and C4(S). Additionally, C4Ivestigator reports the full overall C4A and C4B aligned sequence, enabling nucleotide level analysis. To demonstrate the utility of this workflow we have analyzed C4A and C4B variation in the 1000 Genomes Project Data set, showing that these genes are highly poly-allelic with many variants that have the potential to impact C4 protein function.

High-resolution DNA methylation screening of the major histocompatibility complex in multiple sclerosis.

Background: The HLA-DRB1 gene in the major histocompatibility complex (MHC) region in chromosome 6p21 is the strongest genetic factor identified as influencing multiple sclerosis (MS) susceptibility. DNA methylation changes associated with MS have been consistently detected at the MHC region. However, understanding the full scope of epigenetic regulations of the MHC remains incomplete, due in part to the limited coverage of this region by standard whole genome bisulfite sequencing or array-based methods. Methods: We developed and validated an MHC capture protocol coupled with bisulfite sequencing and conducted a comprehensive analysis of the MHC methylation landscape in blood samples from 147 treatment naïve MS study participants and 129 healthy controls. Results: We identified 132 differentially methylated region (DMRs) within MHC region associated with disease status. The DMRs overlapped with established MS risk loci. Integration of the MHC methylome with human leukocyte antigen (HLA) genetic data indicate that the methylation changes are significantly associated with HLA genotypes. Using DNA methylation quantitative trait loci (mQTL) mapping and the causal inference test (CIT), we identified 643 cis-mQTL-DMRs paired associations, including 71 DMRs possibly mediating causal relationships between 55 single nucleotide polymorphisms (SNPs) and MS risk. Results: The results describe MS-associated methylation changes in MHC region and highlight the association between HLA genotypes and methylation changes. Results from the mQTL and CIT analyses provide evidence linking MHC region variations, methylation changes, and disease risk for MS.

The landscape of the immunoglobulin repertoire in endemic pemphigus foliaceus.

Introduction: Primarily driven by autoreactive B cells, pemphigus foliaceus (PF) is an uncommon autoimmune blistering skin disease of sporadic occurrence worldwide. However, PF reaches a prevalence of 3% in the endemic areas of Brazil, the highest ever registered for any autoimmune disease, which indicates environmental factors influencing the immune response in susceptible individuals. We aimed to provide insights into the immune repertoire of patients with PF living in the endemic region of the disease, compared to healthy individuals from the endemic region and a non-endemic area. Methods: We characterized the B-cell repertoire in i) nontreated patients (n=5); ii) patients under immunosuppressive treatment (n=5); iii) patients in remission without treatment (n=6); and two control groups iv) from the endemic (n=6) and v) non-endemic areas in Brazil (n=4). We used total RNA extracted from peripheral blood mononuclear cells and performed a comprehensive characterization of the variable region of immunoglobulin heavy chain (IGH) in IgG and IgM using next-generation sequencing. Results: Compared to individuals from a different area, we observed remarkably lower clonotype diversity in the B-cell immune repertoire of patients and controls from the endemic area (p < 0.02), suggesting that the immune repertoire in the endemic area is under geographically specific and intense environmental pressure. Moreover, we observed longer CDR3 sequences in patients, and we identified differential disease-specific usage of IGHV segments, including increased IGHV3-30 and decreased IGHV3-23 in patients with active disease (p < 0.04). Finally, our robust network analysis discovered clusters of CDR3 sequences uniquely observed in patients with PF. Discussion: Our results indicate that environmental factors, in addition to disease state, impact the characteristics of the repertoire. Our findings can be applied to further investigation of the environmental factors that trigger pemphigus and expand the knowledge for identifying new targeted and more effective therapies.

High-throughput complement component 4 genomic sequence analysis with C4Investigator.

The complement component 4 gene locus, composed of the C4A and C4B genes and located on chromosome 6, encodes for C4 protein, a key intermediate in the classical and lectin pathways of the complement system. The complement system is an important modulator of immune system activity and is also involved in the clearance of immune complexes and cellular debris. The C4 gene locus exhibits copy number variation, with each composite gene varying between 0-5 copies per haplotype, C4 genes also vary in size depending on the presence of the HERV retrovirus in intron 9, denoted by C4(L) for long-form and C4(S) for short-form, which modulates expression and is found in both C4A and C4B. Additionally, human blood group antigens Rodgers and Chido are located on the C4 protein, with the Rodger epitope generally found on C4A protein, and the Chido epitope generally found on C4B protein. C4 copy number variation has been implicated in numerous autoimmune and pathogenic diseases. Despite the central role of C4 in immune function and regulation, high-throughput genomic sequence analysis of C4 variants has been impeded by the high degree of sequence similarity and complex genetic variation exhibited by these genes. To investigate C4 variation using genomic sequencing data, we have developed a novel bioinformatic pipeline for comprehensive, high-throughput characterization of human C4 sequence from short-read sequencing data, named C4Investigator. Using paired-end targeted or whole genome sequence data as input, C4Investigator determines gene copy number for overall C4, C4A, C4B, C4(Rodger), C4(Ch), C4(L), and C4(S), additionally, C4Ivestigator reports the full overall C4 aligned sequence, enabling nucleotide level analysis of C4. To demonstrate the utility of this workflow we have analyzed C4 variation in the 1000 Genomes Project Dataset, showing that the C4 genes are highly poly-allelic with many variants that have the potential to impact C4 protein function.

A common allele of HLA is associated with asymptomatic SARS-CoV-2 infection.

Studies have demonstrated that at least 20% of individuals infected with SARS-CoV-2 remain asymptomatic1-4. Although most global efforts have focused on severe illness in COVID-19, examining asymptomatic infection provides a unique opportunity to consider early immunological features that promote rapid viral clearance. Here, postulating that variation in the human leukocyte antigen (HLA) loci may underly processes mediating asymptomatic infection, we enrolled 29,947 individuals, for whom high-resolution HLA genotyping data were available, in a smartphone-based study designed to track COVID-19 symptoms and outcomes. Our discovery cohort (n = 1,428) comprised unvaccinated individuals who reported a positive test result for SARS-CoV-2. We tested for association of five HLA loci with disease course and identified a strong association between HLA-B*15:01 and asymptomatic infection, observed in two independent cohorts. Suggesting that this genetic association is due to pre-existing T cell immunity, we show that T cells from pre-pandemic samples from individuals carrying HLA-B*15:01 were reactive to the immunodominant SARS-CoV-2 S-derived peptide NQKLIANQF. The majority of the reactive T cells displayed a memory phenotype, were highly polyfunctional and were cross-reactive to a peptide derived from seasonal coronaviruses. The crystal structure of HLA-B*15:01-peptide complexes demonstrates that the peptides NQKLIANQF and NQKLIANAF (from OC43-CoV and HKU1-CoV) share a similar ability to be stabilized and presented by HLA-B*15:01. Finally, we show that the structural similarity of the peptides underpins T cell cross-reactivity of high-affinity public T cell receptors, providing the molecular basis for HLA-B*15:01-mediated pre-existing immunity.

Comparison between qPCR and RNA-seq reveals challenges of quantifying HLA expression.

Human leukocyte antigen (HLA) class I and II loci are essential elements of innate and acquired immunity. Their functions include antigen presentation to T cells leading to cellular and humoral immune responses, and modulation of NK cells. Their exceptional influence on disease outcome has now been made clear by genome-wide association studies. The exons encoding the peptide-binding groove have been the main focus for determining HLA effects on disease susceptibility/pathogenesis. However, HLA expression levels have also been implicated in disease outcome, adding another dimension to the extreme diversity of HLA that impacts variability in immune responses across individuals. To estimate HLA expression, immunogenetic studies traditionally rely on quantitative PCR (qPCR). Adoption of alternative high-throughput technologies such as RNA-seq has been hampered by technical issues due to the extreme polymorphism at HLA genes. Recently, however, multiple bioinformatic methods have been developed to accurately estimate HLA expression from RNA-seq data. This opens an exciting opportunity to quantify HLA expression in large datasets but also brings questions on whether RNA-seq results are comparable to those by qPCR. In this study, we analyze three classes of expression data for HLA class I genes for a matched set of individuals: (a) RNA-seq, (b) qPCR, and (c) cell surface HLA-C expression. We observed a moderate correlation between expression estimates from qPCR and RNA-seq for HLA-A, -B, and -C (0.2 ≤ rho ≤ 0.53). We discuss technical and biological factors which need to be accounted for when comparing quantifications for different molecular phenotypes or using different techniques.

Genetic Association and Differential RNA Expression of Histone (De)Acetylation-Related Genes in Pemphigus Foliaceus-A Possible Epigenetic Effect in the Autoimmune Response.

Pemphigus foliaceus (PF) is an autoimmune skin blistering disease characterized by antidesmoglein-1 IgG production, with an endemic form (EPF) in Brazil. Genetic and epigenetic factors have been associated with EPF, but its etiology is still not fully understood. To evaluate the genetic association of histone (de)acetylation-related genes with EPF susceptibility, we evaluated 785 polymorphisms from 144 genes, for 227 EPF patients and 194 controls. Carriers of HDAC4_rs4852054*A were more susceptible (OR = 1.79, p = 0.0038), whereas those with GSE1_rs13339618*A (OR = 0.57, p = 0.0011) and homozygotes for PHF21A_rs4756055*A (OR = 0.39, p = 0.0006) were less susceptible to EPF. These variants were not associated with sporadic PF (SPF) in German samples of 75 SPF patients and 150 controls, possibly reflecting differences in SPF and EPF pathophysiology. We further evaluated the expression of histone (de)acetylation-related genes in CD4+ T lymphocytes, using RNAseq. In these cells, we found a higher expression of KAT2B, PHF20, and ZEB2 and lower expression of KAT14 and JAD1 in patients with active EPF without treatment compared to controls from endemic regions. The encoded proteins cause epigenetic modifications related to immune cell differentiation and cell death, possibly affecting the immune response in patients with PF.

Genetic Associations and Differential mRNA Expression Levels of Host Genes Suggest a Viral Trigger for Endemic Pemphigus Foliaceus.

The long search for the environmental trigger of the endemic pemphigus foliaceus (EPF, fogo selvagem) has not yet resulted in any tangible findings. Here, we searched for genetic associations and the differential expression of host genes involved in early viral infections and innate antiviral defense. Genetic variants could alter the structure, expression sites, or levels of the gene products, impacting their functions. By analyzing 3063 variants of 166 candidate genes in 227 EPF patients and 194 controls, we found 12 variants within 11 genes associated with differential susceptibility (p < 0.005) to EPF. The products of genes TRIM5, TPCN2, EIF4E, EIF4E3, NUP37, NUP50, NUP88, TPR, USP15, IRF8, and JAK1 are involved in different mechanisms of viral control, for example, the regulation of viral entry into the host cell or recognition of viral nucleic acids and proteins. Only two of nine variants were also associated in an independent German cohort of sporadic PF (75 patients, 150 controls), aligning with our hypothesis that antiviral host genes play a major role in EPF due to a specific virus−human interaction in the endemic region. Moreover, CCL5, P4HB, and APOBEC3G mRNA levels were increased (p < 0.001) in CD4+ T lymphocytes of EPF patients. Because there is limited or no evidence that these genes are involved in autoimmunity, their crucial role in antiviral responses and the associations that we observed support the hypothesis of a viral trigger for EPF, presumably a still unnoticed flavivirus. This work opens new frontiers in searching for the trigger of EPF, with the potential to advance translational research that aims for disease prevention and treatment.

HLA variation and antigen presentation in COVID-19 and SARS-CoV-2 infection.

The extraordinary variation of the human leukocyte antigen (HLA) molecules is critical for diversifying antigen presentation to T cells. Coupled with the rise of novel strains and rapidly evolving immune evasion by SARS-CoV-2 proteins, HLA-mediated immunity in COVID-19 is critically important but far from being fully understood. A growing number of studies have found the association of HLA variants with different COVID-19 outcomes and that HLA genotypes associate with differential immune responses against SARS-CoV-2. Prediction studies have shown that mutations in multiple viral strains, most concentrated in the Spike protein, affect the affinity between these mutant peptides and HLA molecules. Understanding the impact of this variation on T-cell responses is critical for comprehending the immunogenic mechanisms in both natural immunity and vaccine development.

Population structure and forensic genetic analyses in Guarani and Kaingang Amerindian populations from Brazil.

When DNA profile comparisons between a crime scene trace and a reference sample generate correspondence, the match probability has to be estimated, so that evaluation of the strength of the forensic DNA evidence can be made. The random match probability estimations require information on allele frequencies and an adjustment factor, referred to as theta (θ) or Fst, a co-ancestry correction factor for subpopulation effects. The θ value has been standardized for urban and isolated populations, but inconsistencies have been reported when it is specifically calculated for smaller and isolated populations, including Amerindian populations. Notably, attempts to characterize forensic markers of these minor populations have been extensively limited and more conservative estimates of the correction factor may be generated for each of them. Therefore, we estimate allele frequencies of 21 autosomal STR markers used for forensic testing and calculated relevant forensic parameters for the set. In addition, we featured the possible structure of five Brazilian Amerindian populations that have been genetically isolated for centuries so we could obtain the appropriate θ value for them. The sample consisted of 319 individuals: (1) 121 Kaingang, from two communities: Ivaí (KIV=61) and Rio das Cobras (KRC=60); and (2) 198 Guaranis from three communities: Mbya from Rio das Cobras (GRC=51), Guarani Ñandeva (GND=71) and Guarani Kaiowá (GKW=76). Between Guarani populations low (Rst=0.0402, p < 10-4) to high (Rst=0.1557, p < 10-5) differentiation was found. Regarding Guarani and Kaingang populations, intermediate (Rst=0.0590, p < 10-5) to high (Rst=0.1604, p < 10-5) differentiation was found. The two Kaingang populations showed very low differentiation between them (Rst=0.0017, p = 0.27), which justifies the union of both genetic data for forensic databases and calculations. The combined power of discrimination (PD) and the combined power of exclusion (PE) were calculated for each population, demonstrating the usefulness of this set of markers in forensic and kinship analysis regarding these populations. Considering the demographic heterogeneity of Amerindian populations in general, the Fst mean value (0.03) was evaluated regarding 43 different indigenous populations from the Americas, including Guaranis and Kaingangs. This result confirms the adequacy of the standardized θ value for the forensic random match probability estimations involving Amerindian populations.

Uniparental markers reveal new insights on subcontinental ancestry and sex-biased admixture in Brazil.

The Brazilian population is a product of asymmetric admixture among European men and Amerindian and African women. However, Brazilian subcontinental ancestry is scarcely documented, especially regarding its African roots. Here, we aimed to unveil the uniparental continental and subcontinental contributions from distinct Brazilian regions, including South (n = 43), Southeast (n = 71), the poorly genetically characterized Central-Western region (n = 323), and a subset of unique Brazilian Amerindians (n = 24), in the context of their genome-wide ancestral contributions. The overwhelming majority of European Y haplogroups (85%) contrast sharply with the predominant African and Amerindian mtDNA haplogroups (73.2%) in admixed populations, whereas in Amerindians, non-Native haplogroups could only be detected through the paternal line. Our in-depth investigation of uniparental markers showed signals of an Andean and Central-Brazilian Amerindian maternal contribution to Southeastern and Central-Western Brazil (83.1 ± 2.1% and 56.9 ± 0.2%, respectively), the last having the highest paternal Amerindian ancestry yet described for an admixed Brazilian region (9.7%) and contrasting with higher Southern-Brazilian Amerindian contribution to Southern Brazil (59.6 ± 1%). Unlike the higher African Bantu contribution previously reported for the South and Southeast, a relevant Western African non-Bantu contribution was detected in those regions (85.7 ± 5% and 71.8 ± 10.8% respectively). In contrast, a higher Bantu contribution was described for the first time in the Central-West (64.8 ± 1.3% maternal and 86.9 ± 9.6% paternal). We observed sex-biased signatures consistent with the historically recorded Brazilian colonization and added new insights in the subcontinental maternal ancestry of Brazilians from regions never studied at this level.

Multiple sclerosis therapies differentially affect SARS-CoV-2 vaccine-induced antibody and T cell immunity and function.

BACKGROUNDVaccine-elicited adaptive immunity is a prerequisite for control of SARS-CoV-2 infection. Multiple sclerosis (MS) disease-modifying therapies (DMTs) differentially target humoral and cellular immunity. A comprehensive comparison of the effects of MS DMTs on SARS-CoV-2 vaccine-specific immunity is needed, including quantitative and functional B and T cell responses.METHODSSpike-specific Ab and T cell responses were measured before and following SARS-CoV-2 vaccination in a cohort of 80 study participants, including healthy controls and patients with MS in 6 DMT groups: untreated and treated with glatiramer acetate (GA), dimethyl fumarate (DMF), natalizumab (NTZ), sphingosine-1-phosphate (S1P) receptor modulators, and anti-CD20 mAbs. Anti-spike-Ab responses were assessed by Luminex assay, VirScan, and pseudovirus neutralization. Spike-specific CD4+ and CD8+ T cell responses were characterized by activation-induced marker and cytokine expression and tetramer.RESULTSAnti-spike IgG levels were similar between healthy control participants and patients with untreated MS and those receiving GA, DMF, or NTZ but were reduced in anti-CD20 mAb- and S1P-treated patients. Anti-spike seropositivity in anti-CD20 mAb-treated patients was correlated with CD19+ B cell levels and inversely correlated with cumulative treatment duration. Spike epitope reactivity and pseudovirus neutralization were reduced in anti-CD20 mAb- and S1P-treated patients. Spike-specific CD4+ and CD8+ T cell reactivity remained robust across all groups, except in S1P-treated patients, in whom postvaccine CD4+ T cell responses were attenuated.CONCLUSIONThese findings from a large cohort of patients with MS exposed to a wide spectrum of MS immunotherapies have important implications for treatment-specific COVID-19 clinical guidelines.FUNDINGNIH grants 1K08NS107619, K08NS096117, R01AI159260, R01NS092835, R01AI131624, and R21NS108159; NMSS grants TA-1903-33713 and RG1701-26628; Westridge Foundation; Chan Zuckerberg Biohub; Maisin Foundation.

Inhibition of extracellular traps by spores of Trichoderma stromaticum on neutrophils obtained from human peripheral blood.

Although the genus Trichoderma is widely used as a biocontrol agent in crops, little is known about its potential impact on the human immune system. In mice, our group has shown that exposition to T. asperelloides spores lead to reduced neutrophil counts in the peripheral blood and in the peritoneal cavity. In addition, T. stromaticum spores produced an inflammatory infiltrate on mice lungs, reducing the levels of IFN-γ and IL-10 cytokines, reactive oxygen species, and receptors of microbial patterns. Here we demonstrate that the interaction of human peripheral neutrophils with T. stromaticum spores also leads to a reduced release of neutrophil extracellular traps (NETs) after induction with the NET-inducer agent phorbol 12-myristate 13-acetate. This interaction also reduced the expression levels of multiple microRNAs, such as miR-221, miR-222, miR-223 and miR-27a, as well as genes related to NETs, such as ELANE, MPO and PADI4. Furthermore, T. stromaticum spores affected the expression of the genes SOCS3, TLR4, CSNK2A1, GSDMD, and NFFKBIA, related to the activation of inflammatory immune responses in neutrophils. Overall, our results suggest T. stromaticum as a potential NET inhibitor and as an immunomodulatory agent. Since this fungus is used as biocontrol in crops, our findings point to the importance of advancing our knowledge on the effects of this bioagent on the human immune system. Finally, the study of the active compounds produced by the fungus is also important for the prospection of new drugs that could be used to block the exacerbation of inflammatory immune responses present in several human diseases.

A common allele of HLA mediates asymptomatic SARS-CoV-2 infection.

Despite some inconsistent reporting of symptoms, studies have demonstrated that at least 20% of individuals infected with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) will remain asymptomatic. Although most global efforts have focused on understanding factors underlying severe illness in COVID-19 (coronavirus disease of 2019), the examination of asymptomatic infection provides a unique opportunity to consider early disease and immunologic features promoting rapid viral clearance. Owing to its critical role in the immune response, we postulated that variation in the human leukocyte antigen (HLA) loci may underly processes mediating asymptomatic infection. We enrolled 29,947 individuals registered in the National Marrow Donor Program for whom high-resolution HLA genotyping data were available in the UCSF Citizen Science smartphone-based study designed to track COVID-19 symptoms and outcomes. Our discovery cohort (n=1428) was comprised of unvaccinated, self-identified subjects who reported a positive test result for SARS-CoV-2. We tested for association of five HLA loci (HLA-A, -B, -C, -DRB1, -DQB1) with disease course and identified a strong association of HLA-B*15:01 with asymptomatic infection, and reproduced this association in two independent cohorts. Suggesting that this genetic association is due to pre-existing T-cell immunity, we show that T cells from pre-pandemic individuals carrying HLA-B*15:01 were reactive to the immunodominant SARS-CoV-2 S-derived peptide NQKLIANQF, and 100% of the reactive cells displayed memory phenotype. Finally, we characterize the protein structure of HLA-B*15:01-peptide complexes, demonstrating that the NQKLIANQF peptide from SARS-CoV-2, and the highly homologous NQKLIANAF from seasonal coronaviruses OC43-CoV and HKU1-CoV, share similar ability to be stabilized and presented by HLA-B*15:01, providing the molecular basis for T-cell cross-reactivity and HLA-B*15:01-mediated pre-existing immunity.

Remarkably Low KIR and HLA Diversity in Amerindians Reveals Signatures of Strong Purifying Selection Shaping the Centromeric KIR Region.

The killer-cell immunoglobulin-like receptors (KIR) recognize human leukocyte antigen (HLA) molecules to regulate the cytotoxic and inflammatory responses of natural killer cells. KIR genes are encoded by a rapidly evolving gene family on chromosome 19 and present an unusual variation of presence and absence of genes and high allelic diversity. Although many studies have associated KIR polymorphism with susceptibility to several diseases over the last decades, the high-resolution allele-level haplotypes have only recently started to be described in populations. Here, we use a highly innovative custom next-generation sequencing method that provides a state-of-art characterization of KIR and HLA diversity in 706 individuals from eight unique South American populations: five Amerindian populations from Brazil (three Guarani and two Kaingang); one Amerindian population from Paraguay (Aché); and two urban populations from Southern Brazil (European and Japanese descendants from Curitiba). For the first time, we describe complete high-resolution KIR haplotypes in South American populations, exploring copy number, linkage disequilibrium, and KIR-HLA interactions. We show that all Amerindians analyzed to date exhibit the lowest numbers of KIR-HLA interactions among all described worldwide populations, and that 83-97% of their KIR-HLA interactions rely on a few HLA-C molecules. Using multiple approaches, we found signatures of strong purifying selection on the KIR centromeric region, which codes for the strongest NK cell educator receptors, possibly driven by the limited HLA diversity in these populations. Our study expands the current knowledge of KIR genetic diversity in populations to understand KIR-HLA coevolution and its impact on human health and survival.

Population-specific diversity of the immunoglobulin constant heavy G chain (IGHG) genes.

Human immunoglobulin G (IgG) molecules, IgG1, IgG2 and IgG3, exhibit substantial inter-individual variation in their constant heavy chain regions, as discovered by serological methods. This polymorphism is encoded by the IGHG1, IGHG2, and IGHG3 genes and may influence antibody function. We sequenced the coding fragments of these genes in 95 European Americans, 94 African Americans, and 94 Black South Africans. Striking differences were observed between the population groups, including extremely low amino acid sequence variation in IGHG1 among South Africans, and higher IGHG2 and IGHG3 diversity in individuals of African descent compared to individuals of European descent. Molecular definition of the loci illustrates a greater level of allelic polymorphism than previously described, including the presence of common IGHG2 and IGHG3 variants that were indistinguishable serologically. Comparison of our data with the 1000 Genome Project sequences indicates overall agreement between the datasets, although some inaccuracies in the 1000 Genomes Project are likely. These data represent the most comprehensive analysis of IGHG polymorphisms across major populations, which can now be applied to deciphering their functional impact.

Impact of multiple sclerosis disease-modifying therapies on SARS-CoV-2 vaccine-induced antibody and T cell immunity.

Vaccine-elicited adaptive immunity is an essential prerequisite for effective prevention and control of coronavirus 19 (COVID-19). Treatment of multiple sclerosis (MS) involves a diverse array of disease-modifying therapies (DMTs) that target antibody and cell-mediated immunity, yet a comprehensive understanding of how MS DMTs impact SARS-CoV-2 vaccine responses is lacking. We completed a detailed analysis of SARS-CoV-2 vaccine-elicited spike antigen-specific IgG and T cell responses in a cohort of healthy controls and MS participants in six different treatment categories. Two specific DMT types, sphingosine-1-phosphate (S1P) receptor modulators and anti-CD20 monoclonal antibodies (mAb), resulted in significantly reduced spike-specific IgG responses. Longer duration of anti-CD20 mAb treatment prior to SARS-CoV-2 vaccination were associated with absent antibody responses. Except for reduced CD4+ T cell responses in S1P-treated patients, spike-specific CD4+ and CD8+ T cell reactivity remained robust across all MS treatment types. These findings have important implications for clinical practice guidelines and vaccination recommendations in MS patients and other immunosuppressed populations.

An immunogenetic view of COVID-19.

Meeting the challenges brought by the COVID-19 pandemic requires an interdisciplinary approach. In this context, integrating knowledge of immune function with an understanding of how genetic variation influences the nature of immunity is a key challenge. Immunogenetics can help explain the heterogeneity of susceptibility and protection to the viral infection and disease progression. Here, we review the knowledge developed so far, discussing fundamental genes for triggering the innate and adaptive immune responses associated with a viral infection, especially with the SARS-CoV-2 mechanisms. We emphasize the role of the HLA and KIR genes, discussing what has been uncovered about their role in COVID-19 and addressing methodological challenges of studying these genes. Finally, we comment on questions that arise when studying admixed populations, highlighting the case of Brazil. We argue that the interplay between immunology and an understanding of genetic associations can provide an important contribution to our knowledge of COVID-19.

High-throughput Interpretation of Killer-cell Immunoglobulin-like Receptor Short-read Sequencing Data with PING.

The killer-cell immunoglobulin-like receptor (KIR) complex on chromosome 19 encodes receptors that modulate the activity of natural killer cells, and variation in these genes has been linked to infectious and autoimmune disease, as well as having bearing on pregnancy and transplant outcomes. The medical relevance and high variability of KIR genes makes short-read sequencing an attractive technology for interrogating the region, providing a high-throughput, high-fidelity sequencing method that is cost-effective. However, because this gene complex is characterized by extensive nucleotide polymorphism, structural variation including gene fusions and deletions, and a high level of homology between genes, its interrogation at high resolution has been thwarted by bioinformatic challenges, with most studies limited to examining presence or absence of specific genes. Here, we present the PING (Pushing Immunogenetics to the Next Generation) pipeline, which incorporates empirical data, novel alignment strategies and a custom alignment processing workflow to enable high-throughput KIR sequence analysis from short-read data. PING provides KIR gene copy number classification functionality for all KIR genes through use of a comprehensive alignment reference. The gene copy number determined per individual enables an innovative genotype determination workflow using genotype-matched references. Together, these methods address the challenges imposed by the structural complexity and overall homology of the KIR complex. To determine copy number and genotype determination accuracy, we applied PING to European and African validation cohorts and a synthetic dataset. PING demonstrated exceptional copy number determination performance across all datasets and robust genotype determination performance. Finally, an investigation into discordant genotypes for the synthetic dataset provides insight into misaligned reads, advancing our understanding in interpretation of short-read sequencing data in complex genomic regions. PING promises to support a new era of studies of KIR polymorphism, delivering high-resolution KIR genotypes that are highly accurate, enabling high-quality, high-throughput KIR genotyping for disease and population studies.